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ABSTRACT 

The distribution of i^C^^O J = 9 ^ 8 (1.037 THz) emission has been mapped in OMC-1 at 35 points 
with 84" resolution. This is the first map of this source in this transition and only the second velocity- 
resolved ground-based observation of a line in the THz frequency band. There is emission present at 
all points in the map, a region roughly 4' x 6' in size, with peak antenna temperature dropping only 
near the edges. Away from the Orion KL outflow, the velocity structure suggests that most of the 
emission comes from the OMC-1 photon-dominated region, with a typical line width of 3-6 km s~^. 
Large velocity gradient modeling of the emission in J = 9 — > 8 and six lower transitions suggests that 
the lines originate in regions with temperatures around 120 K and densities of at least 10^'^ cm~'^ near 
O^C Ori and at the Orion bar, and from 70 K gas at around 10^ cm~^ southeast and west of the bar. 
These observations are among the first made with the 0.8 m Smithsonian Astrophysical Observatory 
Receiver Lab Telescope, a new instrument designed to observe at frequencies above 1 THz from an 
extremely high and dry site in northern Chile. 

Subject headings: ISM: clouds - ISM: individual: OMC-1 - ISM: individual: Orion Klcinmann-Low - 
ISM: individual: Orion Bar - submillimeter - telescopes 



1. introduction 

The OMC-1 region of the Orion A molecular cloud 
is the nearest site (^^500 pc) of recent high mass star 
formation. Its diverse components make it a natural 
laboratory for the study of many stages of star forma- 
tion. The OMC-1 cloud is one of three clouds connected 
by the Orion "ridge." a l ong and dense molecular fila- 
ment ijCastets et al.lll990|) . The H II regions M42 and 
M43 sit in front of the ridge, and extended narrow-line 
CO emission (the "spike" component) is excited in the 
photon-dominated region (PDR) at the interface between 
the ionized and neutral gas. The brightest star in the 
Trapezium, 6^G Ori, is the dominant source for the far- 
ultraviolet photons t hat sustain the M4 2 H II region and 
its associated PDR (jStacev et al.lll99^ . One arcminute 
to the north-west oi O^C Ori, the Kleinmann-Low nebula 
(KL) is the brightest source outside the solar svstem at 
20 /im, with a luminosity of ~ lO^L© (jKleinmann fc Low! 
IT967i^ . Several powerful sources have been identi- 
fied in the infrared within KL, including the Becklin- 
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Neug ebauer point source (BN; IBecklin fc Neugebaueil 
1967l i. likely a runaway B star ijPlambeck et alJ Il995t 



TanI l2003) . and IRc2 once thought to be the power 
source for the nebula (|Downes et al.llT98ll) but more re- 
cently reso lved into several sourc e s that may not be self- 
luminous (iPougados et"an 119931 iMenten fc ReidI 119951 
iCezari et al.lll998h . Molecular line observations of Orion 
KL show emission extending out to nearly ±100 km s~^ 
(e.g. CO J = 7 ^ 6 and J = 4 ^ 3 in iWilson et all 
I2001L hereafter WMKHOl). This large line width, usu- 
ally called "platea u" emission, is due to an outflow long 
attributed to IRc2 ijWright et alJl98^ that is more likely 
associated with nearb y radio and infr ared p oint sources 
(Menten fc Reid 19951: iGreen hiU et aT] ll99ah . The Orion 
bar is a well studied region of the PDR south-east of the 
Trapezium. At this location the ionization front is viewed 
edge-on over a length of ~3', allowing direct examina- 
tion of the chemical s tratification in a plan e-parallel P DR 
(e-g.- lvan der Werf e t al. 1996; Walmslev etani200l . 

The OMC-1 cloud has been mapped in many 
millime ter, submill i meter, and far-IR transitions 
(e.g.. iStacev et alJ Il993t IC oldsmith et al. 199^ 
i Herrniann et al.lll997t lPlume"et'al. . .200Q: .Sempere et all 



I. Previ ous observa- 
tions in CO transitions up to J = 7 — > 6 l)Schiilz et alJ 
JJ.95; WMKHOl) have found warm PDR gas extending 
over a very large region. In this paper we present 
the first map of OMC-1 in the J = 9 — > 8 transition 
of ^^C^^O. This line has previously been obser ved in 
OMC-1 from the Kuiper Airborne Observatory (fRoserl 
11991; Betz & Boreiko 1993) and fr om the Heinrich Hertz 
Submillimeter Telescope l(Kawamtira et all I2002D . but 
never at more than a few points. Our map contains 35 
pointings and covers ~25 arcmin^ at 84" resolution. Our 
spectra provide information about the population of the 
J = 9 rotational level (250 K), nearly 100 K higher than 
the J = 7 level probed by J = 7 ^ 6 maps, and the 
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critical density at J 9 ^ 8 10® cm^"^) is approxi- 
mately twice that atJ = 7^6(^5x 10^ cm^'^). This 
THz transition should therefore be a better probe of 
the density and temperature of warm and dense sources 
such as OMC-1, where lower transitions are thermalized. 

Because of absorption by the Earth's atmosphere, the 
J = 9 ^ 8 transition of CO as well as other astro- 
nomically interesting lines at frequencies above 1 THz 
are unobservable from current ground-based observato- 
ries except in very unusual weather. However, recent at- 
mospheric measurements have demonstrated that several 
windows can open up between 1 and 3 THz at very high 
and dry locations (jPaine et al.l l200fi iMatsiishita et all 
^99). Under favorable conditions, three of these win- 
dows, centered at 1.03, 1.35, and 1.5 THz, show trans- 
mission as high as 40%. These windows include six ro- 
tational transitions of ^^CO and "'^^CO, as well as the 
ground-state fine-structure transition of singly ionized 
nitrogen, one of the most important coolants of the dif- 
fuse ISM, and numerous other common and exotic atomic 
and molecular lines. A new telescope, the Receiver Lab 
Telescope (RLT), has been constructed by the Smith- 
sonian Astrophysical Observatory to take advantage of 
these windows, and is currently deployed at 5525 m alti- 
tude in northern Chile. This map of OMC-1 represents 
some of the first observations made with the RLT and 
demonstrates the capabilities of this new instrument. 

2. OBSERVATIONS 

We observed the OMC-1 region of the Orion A molec- 
ular cloud in the J = 9 ^ 8 transition of ^^C^^O using 
the 0.8 m diameter RLT (glundeU et al. 2002) on Cerro 
Sairecabur, Chile. At the 1036.912 GHz frequency of this 
line we expect a beamsize of 84". Our 330 channel auto- 
correlation spectrometer has a bandwidth of 286 km s~^ 
(990 MHz) and a resolution of 1.04 km s^^ at this fre- 
quency. The data were acquired on the morning of 2002 
December 15 between 05:30 and 07:30 UT. We obtained 
spectra at 35 points with 50" spacing in an irregularly 
shaped region of OMC-1. These observations were made 
in position-switching mode with a reference point 10' 
to the west of each map point. Each spectrum is the 
sum of 12 integration cycles of 5 s on-source and 5 s 
off-source. The map is centered on a = 05''35'"14'' and 
S = -05°22'31" (J2000.0), which coincides with the lo- 
cation of several Orio n line surv eys (Sutton et al. 1985; 
iBlake et aLl lTgSe: Schil ke et all 12001(1 . Orion KL, BN, 
and IRc2 are all well within the beam at this location. 
The ranges in offsets from the map center are Aa = —50" 
to +150" and A<5 = -250" to +100". 

The telescope pointing varies from night to night by 
up to tens of arcseconds because of shifts in the tele- 
scope foundation. To remove these variations we used 
an optical guidescope to determine offsets from our base 
pointing model. The limitation of this method is the 
accuracy with which we know the offset between the op- 
tical and radio beams, about 15", and we take this value 
to be our absolute pointing error. Within the map the 
relative pointing is considerably better. To ensure that 
there was no systematic drift in the (0", 0") point over 
the observation period we remeasured the pointing off- 
sets on a nearby star approximately every 25 minutes. 
We also measured the center spectrum twice; the two 
measurements were separated by 1 hr, and there were 



three measurements of the pointing offsets made during 
that interval. As a check of the consistency of these two 
spectra, we rescale them so that the averages of their 
three highest channels are the same and compute the 
rms difference between them in the nonbaseline channels. 
When the spectra are scaled to 100 K we measure an rms 
difference of 4.8 K, very similar to the 4.2 K we expect 
from the noise seen in the baseline channels, and the dif- 
ference spectrum shows no remnant structure. Adjacent 
spectra are significantly different; the most similar spec- 
trum, at (0", —50"), differs from the two center-point 
spectra by 7.7 and 7.4 K and has narrower wings. We 
also note that the averages of the three highest channels 
are 65.5 and 64.1 K for the first and second center-point 
spectra, respectively, indicating that the calibration was 
self-consistent to 2% over this interval. 

The hot electron bolometer recei ver used for these mea - 
surements is described elsewhere l|Meledin et al.ll2004|) . 
During these observations the double-sideband receiver 
temperature was found to be 980 K through observations 
of cold (liquid nitrogen) and ambient loads. Atmospheric 
transmission measurements were mad e every 10 minutes 
by a Fourier transform spectrometer ijPaine et al.ll20d(Tl 
located a few meters away from the telescope. Over 
the course of the observations the zenith transmission 
at the line frequency decreased from 23.5% to 18%. To 
determine a single-sideband system temperature (Tsys), 
we assume equal gain in the signal and image sidebands 
(centered on 1036.912 and 1030.848 GHz, respectively) 
because the mixer frequency response is essentially flat 
between 900 GHz and 1.05 THz. Our data are cali- 
brated to the T^ scale of lUhch fc Haa s (1976). W e have 
measured a beam efficiencv (rji in lUlich fc Haaslll97^ 
of 45% at 883 GHz by mapping a telluric ozone ab- 
sorption feature against Jupiter (angular diameter 42"; 
^arronc ct al. 2004). This efficiency is much lower than 
expected; the RLT primary mirror has an rms figure error 
of only 3 fim, causing a 2% Ruze loss, and our blockage 
and illumination pattern give a theoretical efficiency of 
94%. We believe that this indicates that the receiver was 
poorly aligned with the telescope and that a significant 
portion of the receiver beam was being truncated be- 
fore reaching the secondary mirror. We have not made 
this beam measurement nearer to the CO J = 9 ^ 8 
frequency, so we adopt this efficiency for the calibra- 
tion. After accounting for the above factors, the single- 
sideband system temperature at the source elevation in- 
creased from 2.7x lO'' to 6.4x 10^ K during our observa- 
tions as the transmission degraded and Orion set. From 
estimates of the uncertainties in the atmospheric trans- 
mission and the efficiency, we expect that our calibration 
is accurate to ±25%. 

3. RESULTS 

The calibrated spectrum at each of the 35 map points 
is shown in Figure ^ superposed on a contour map of 
the CO J = 9 ^ 8 line fiux integrated between —30 
and -|-40 km s^"'^. The narrow spike emission from the 
PDR is evident throughout the map, while broad plateau 
emission is only observed near the center, around Orion 
KL. Linear baselines were fitted to 27 channels (—47 to 
—30 km s"^ and +40 to -1-47 km s~^) in each spectrum 
and removed. Because of mismatched baselines in the 
output of the three independent sections of our auto- 
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Fig. 1.— CO J = 9 ^ 8 emission in OMC-1. The contours 
trace the hne flux between —30 and +40 km (in K km s~^). 
The contours are spaced by 10% of the peak flux (1260 K km s~^ 
near the center point) and run from 10 to 90%, with the 50% 
contour reinforced. The small crosses mark the location of map 
pointings, and the bold cross marks the center point of the map 
at a = 05'"35'"14=, 5 = -05°22'31" (J2000.0). The individual 
spectra lie on a square grid with 50" spacing. The star marks the 
position of 9^C Ori, and the triangle marks the position of the IRc2 
complex of sources. The numbers 1-4 correspond to the positions 
in Tables HI and 1^ The spectrum at each map point is also plotted. 
The scale for each spectrum is the same as that given in the lower 
left. 



correlation spectrometer, it was only possible to use the 
middle 110 channels (covering —47 to +47 km s~^) in 
our final spectra. Unfortu nately, other observations of 
Orion KL in this transition (IRoseillT991l: IBetz &: Boreikd 
[1993: Kaw amura ct al. 2002) and lower transitions (e.g., 
ISchulz et a l. 1995; WMKHOf) show plateau emission at 
velocities outside the —30 to +40 km s~^ range. Hence, 
the baselines subtracted from the (Aa = 0", = 0") and 
(0", —50") spectra, where the emission is widest in veloc- 
ity, are likely to be inaccurate, and the integrated emis- 
sion obtained from these spectra is significantly dimin- 
ished. 

Figure^shows that the peak emission in the J = 9 ^ 8 
line is ~20" south of our map center. Observations of 
lower CO transition s, from J = 1 ^ and J = 2 ^ 1 
jCastets e t 111119901) to J = 7 ^ 6 l|Howe et al.lll993t 
IWMKHOir show peak emission at the location of IRc2, 
marked with the triangle in Figure The addition of 
reasonable flux corrections at the two central points af- 
fected by the baseline problem can move the center point 
toward (0", 0") by as much as 10". The remaining error 
is consistent with our absolute pointing uncertainty of 
15". 



The outflow near IRc2 is unresolved in our map, as 
expected from its -50" size in J = 7 -> 6 (WMKHOl). 
We have used the emission in the line wings, which is 
spatially confined to this outflow, as a source to map 
the beam and estimate the beam size. For this estimate 
we assume Gaussian profiles for the source and beam, 
the simplest possible assumption. We obtain a range 
of half-power widths between 85" and 105", which, after 
deconvolution with a 50" source, give beam sizes between 
69" and 92". These values provide confirmation that the 
telescope beam is similar to the 84" we predict from the 
telescope and receiver optics design. 

4. DISCUSSION 

4.1. Radiative Transfer Modeling 

CO J = 9 — > 8 emission is present over the entire 
mapped region, and the peak lin e temperature decreases 
only at the edges. iWMKHOll also found J = 7 ^ 6 
emission over their entire (slightly smaller) mapping re- 
gion. As mentioned previously, we expect the J = 9 — > 8 
line to be a better probe of the warm and dense OMC-1 
PDR. To examine the conditions in the PDR, we com- 
pare the emission in this line to that in lower transi- 
tions: i^CO J = 7^6andJ = 4^3 (WMKHOl), 
J = 6 - > 5 (K . N. AUers 2005, in preparation), J = 3 ^ 2 
llTiggeslll99^ ^^CO J = 3 -> 2 (Tiggcs 1993.) . and 
J = 1 ^ (jPlume et alj r200G). The data sets have 
been convolved to the 84" resolution of the RLT, ex- 
cept for the J — 6 ^ 5 data, which have 86" resolution. 
We estimate the density and temperature at various map 
points usi ng the large velocity gradient (LVG) model of 
iStutzki fc Winncwisscr ( 1985). LVG models assume that 
ordered line-of-sight motions are large compared to local 
random velocities and that there is a one-to-one corre- 
spondence between position along the line of sight and 
veloci ty, allowing the radiative transfer to be treated lo- 
cally l)Scoville fc Solomonlll974l) . The free parameters in 
these models are: the molecular hydrogen density [n(H2), 
or simply n], the CO column density per velocity inter- 
val (dN/dv), and the kinetic temperature (Tfci„). The 
model performs calculations t hrough J = 20 with th e 
H2-CO collision coefficients of iFlower fc Launavl l)1985(l . 
Four representative points were selected from the map 
(marked in Figure ^1; the Rayleigh- Jeans main-beam 
brightness temperatures in each transition at each po- 
sition are given in Table ^ To correct the J = 9 ^ 8 
data from T^ to Tmb, we use a source coupling efficiency 
of 85% because of the extended nature of the spike emis- 
sion. This is the same coupling efficiency used with the 
J = 6 ^ 5 data by K. N. AUers (2005, in preparation) for 
a similar beam size. Coupling efficiencies for the other 
data are obtained from the original papers. Position 1 is 
near 9^G Ori, position 2 lies on the Bar, and positions 3 
and 4 fall on the edges of the mapped region where the 
J = 9 ^ 8 emission is decreasing. 

We estimate the model parameters using a mini- 
mization analysis. The peak brightness temperature in 
each line is determined from a Gaussian fit. These tem- 
peratures are compared to a grid of LVG models cover- 
ing Tktn from 40 to 200 K, n from 10^-^ to 10^-^ cm'^, 
and dN/dv from 10^^'^ to 10^° cm~^ s km~^. We as- 
sign an error of 25% to each temperature to account 
for the uncertainty in the inter-telescope calibration and 
the error-beam contribution to the observed temper- 
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atures on this extended source. For calculations of 
^•^CO line temperatures we use a column density ratio 
of ■/V(^^CO)/jV("C O) = 67, as determined in Orion A 
bv lLanger fc Pen'zTa s (1990), and assume that the ^^CO 
and ^'^CO emission arise in gas of the same temperature. 
The Tkin values we derive are very simi lar to those mea- 
sured for ^^CO bv iPlume et al.1 1)2000(1 . where the gas 
kinetic temperature was determined from measurements 
of "'^^CO lines only, suggesting that this assumption is a 
good approximation for this source. 

The parameters derived from these fits are given in Ta- 
ble [3 The goodness of fit at each position is judged by 
examining the minimum value of the reduced (Xr); 
or divided by the number of degrees of freedom. The 
ranges in the table describe the region in parameter space 
where Xr is less than 1 plus the minimum value. The 
observed minima in Xr indicate good fits at all four po- 
sitions. 

The fits at positions 1 and 2 provide examples of the 
difficulties of determining n from sets of thermalized 
transitions. When the ^^CO line temperatures at these 
positions are converted from the Rayleigh- Jeans scale to 
the Planck scale (under the assumption that they are op- 
tically thick), they are roughly constant, indicating LTE 
emission. In this equilibrium case, the line temperatures 
are no longer sensitive to density and the models can 
only set a lower limit to the density of the gas, that be- 
ing the density required to roughly thermalize all of the 
observed transitions. It is only because the temperatures 
of the two ^'^CO transitions are simultaneously matched 
only in a small region of densities that these positions 
have density upper bounds. Higher densities at the same 
column density per velocity interval are not strongly ex- 
cluded because they give only slightly larger values of 
xf- Positions 3 and 4 indicate the value of high- J data 
in dense environments. At these positions the J = 9 — > 8 
emission is not thermalized, allowing strong constraints 
on the density. When we apply our fitting technique to 
these positions without the J = 9 ^ 8 information we 
do not obtain density upper bounds. 

These results match well with previous measurements, 
although we are able to place tighter constraints on the 
gas conditi ons. Based on ^^C O J = 7 ^ 6 and ^^CO 
J = 2 ^ 1, iHowe et al.l l(1993j) found that the CO emis- 
sion could be produced by column densities per veloc- 
ity interval in the range 10^^'^ — lO^*'^ cm~^ s km~^, 
densiti es of 10'^ — 10^ cm~ ^, and temperatures of 40 — 
200 K. iPhime et all ^2m(t) also used LVG models and 
an assumed density of 10 cm~'^ to find a nearly con- 
stant ^'^CO column density per velocity interval of 

JQ16.4±0.3 ^^-^2 g i^j^-l^ ;^ol8.1±0.3 ^^-2 g ^^-1 foj. 

^^CO with their assumed abundance ratio, and temper- 
atures of ~100 K. At all four positions we find a CO 
column within the range 10^^'^ — 10^^'^ cm~^ s km~^ 
and molecular hydrogen density lO'^'^ — 10^'^ cm~^, al- 
though higher densities are allowed at lower likelihood 
at positions 1 and 2. There are significant temperature 
variations over the OMC-1 cloud, with the lowest tem- 
peratures occurring at the two positions farthest from 
the strong UV source 9^C Ori. 

4.2. Line Features 

Away from the broad emission component around 
Orion KL, the quiescent (spike) emission dominates. Fig- 
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Fig. 2.— Velocity of CO J = 9 ^ 8 emission in OMC-1. The line- 
center velocities are obtained from each spectrum using a Gaussian 
fit. A single component is used away from the KL outflow, while 
at the nine center positions a broad component is added but con- 
strained to have the same velocity. Contours trace velocities from 
8 to 10 km in steps of 0.5 km s~^. Symbols are the same as 
in Figure m 



ure|21is a map of the line-center velocities obtained from 
Gaussian fits to the J = 9 ^ 8 spectra. At the nine 
positions closest to the map center where the KL out- 
flow is visible, we add a broad component to the fit but 
constrain its velocity to match that of the narrow compo- 
nent. As seen in lower transitions, such as the J = 7^6 , 
J = 4 ^ 3, and J = 2 ^ 1 transitions in IWMKHOU 
there is a gradient in the line-center velocity of the spike 
emission. The minimum (least redshifted) velocity is ap- 
proximately 7.5 km s~^ near (—50", —100"), increasing 
toward the north and southeast to more than 10 km s~^ 
at the eastern edge of the map. The observed velocities 
are systematically redshifted (typically 0.5—1 km s~^) 
from veloci ties measured in the dense ridg e gas (e.g., CN 
emission in iRodriguez- Franco et aLll200lD . We interpret 
this as evidence that most of the CO J = 9 — > 8 emission 
arises in the PDR rather than gas deeper in the molecu- 
lar cloud. The rcdshift between the PDR surface and the 
ridge may be due to the slow expansion of the H II region 
into the molecular cloud. Emission from the ridge, if it 
is present, could be visible as a blue shoulder in individ- 
ual spectra, but the only spectrum that may show this 
second component is (0", —200"). 

The FWHM of the spike emission, as measured from 
the Gaussian fits, narrows away from the map center, de- 
creasing from 6 km (FWHM) near the KL outflow to 
about 3 km s"""^ in the southeast corner and 2.5 km s""'^ in 
the northeast. WMKHOl shows the line width increasing 
to the west at the declination of 6^C Ori in J = 7 — > 6, 
but our angular resolution does not allow us to disentan- 
gle the KL outflow and compare directly. 
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Only one spectrum shows a clear emission component 
other than the spike and plateau at our resolution and 
signal-to-noise ratio. At the Orion bar [(100", -200"); 
position 2 above], a blue component stands out above 
3 km s^^. The same component may also be present at 
(100", -150") and (50", -200"). The emission at this 
location appears to be the sum of two components: one 
around 10.5 km (consistent with the velocity at the 
Bar), and another around 7.5 km s~^; both components 
are visible in the other ^^CO and ^^CO transitions as 
well. The i^CO and "CO J = 1 ^ channel maps of 
iTauber et al. ( 1994j) resolve the two components in space 
and velocity, with the 7.5 km s~^ component centered 
approximately 30" north and west of the bar emission. 
When convolved to the same angular resolution (84"), 
the line shape is very similar in the J = 4 ^ 3 through 
J = 9 ^ 8 transitions, suggesting that the line-center 
and wing components are at similar average density and 
temperature and differ only in beam filling factor. We 
suggest that the blue component may represent molec- 
ular gas left over from the disrupted near edge of the 
molecular cloud. 

5. CONCLUSIONS 

We have mapped CO J = 9 8 emission in the OMC- 
1 region of the Orion A molecular cloud over a 25 arcmin^ 
region at 84" resolution. The line is detected over the en- 
tire map, similar to what has been seen in other submil- 
limeter CO transitions. The high critical density of this 
transition and the high energy of the J = 9 level make 
it more difficult to excite than lower levels, allowing us 



to measure the molecular hydrogen density, temperature, 
and CO column density per velocity interval of the gas 
at four locations in the PDR by combining these spectra 
with lower transitions. We find that near the center of 
the map even this transition may not be high enough to 
strongly constrain the density because of the high gas 
temperatures. A wider survey of the cloud in this line or 
the higher rotational lines of ^^CO and ^^CO available to 
the RLT should reveal more cool gas for which the THz 
CO transitions can tightly constrain the conditions. 

This map demonstrates that it is possible to make re- 
liable observations from the ground in a region of the 
electromagnetic spectrum normally considered to be in- 
accessible. In the near future we hope to begin obser- 
vations at higher frequencies to make use of the other 
atmospheric windows available from our site. 
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TABLE 1 

Line peak temperatures used in LVG model fitting 



Position^ 


rnl2 b 


rpl2 C 


rpl2 d 


rrl2 C 


rnl2 o 
^3^2 


q^l3 c 
'-3^2 


rpl3 f 


1 


85±21 


137±34 


83±21 


102±26 


93±23 


24±6 


11±3 


2 


80±20 


120±30 




102±26 


86±22 


24±6 


11±3 


3 


36±9 




46±12 




71±18 


12±3 


7.0±1.8 


4 


27±7 




56±14 




74±19 


30±7 


13±3 



Note. — Line peak temperatures labeled as Ty , where X is the carbon isotope 
and Y is the rotational transition. All values are in Kelvins. 
''See labeled positions in Figure [Tl 
''Data from this work. 
<=Data from lWMKHOl 

''Data from K. N. AUers 2005, in preparation. 

<=Data from'figgcs (1993)^ 

'Data from iPlume et al.1 <2000h . 



TABLE 2 

Gas parameters derived from LVG model 



Position Tfei„ log n \og{dN/dv) Xr.min 

(K) (cm^^) (cm^^ s km^') 

1 120±20 3.5 -4.6=^ 18.0 ± 0.2 0.83 

2 120±20 3.5 -4.3=^ 18.0 ± 0.2 0.81 

3 70±10 4.1 ±0.3 17.7 ±0.2 1.46 

4 70±10 3.8 ±0.3 18.1 ±0.2 1.16 



''Poorly constrained upper limit. See text for details. 



